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Abstract 

We give an exposition of recent results on regularity and Fredholm properties for 
first-order one-dimensional hyperbolic PDEs. We show that large classes of bound- 
ary operators cause an effect that smoothness increases with time. This property is 
the key in finding regularizers (parametrices) for hyperbolic problems. We construct 
regularizers for periodic problems for dissipative first-order linear hyperbolic PDEs 
and show that they are modeled by Fredholm operators of index zero. 
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1 Introduction 



In contrast to ODEs and parabolic PDEs, the Fredholm property and regularity behavior 
of hyperbolic problems are much less understood. In a recent series of papers [17, 19, 20], 
the latter two written jointly with Lutz Recke, we undertook a detailed analysis of this 
subject for first-order one-dimensional hyperbolic operators. The purpose of the present 
survey paper is to present some of our results and their extensions with emphasize on the 
smoothing phenomenon, construction of regularizers, and the Fredholmness of index zero. 

An important step in local investigations of nonlinear differential equations (many 
ODEs and parabolic PDEs) is to establish the Fredholm solvability of their linearized 
versions. In the hyperbolic case this step is much more involved. Since the singularities 
of (semi-) linear hyperbolic equations propagate along characteristic curves, a solution 
cannot be more regular in the entire time-space domain than it is on the boundary. 
It can even be less regular which is known as the loss-of-smoothness effect. Therefore 
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the Fredholm analysis of hyperbolic problems requires establishing an optimal regularity 
relation between the spaces of solutions and right-hand sides of the differential equations. 

Proving a Fredholm solvability is typically based on the basic fact that any Fredholm 
operator is exactly a compact perturbation of a bijective operator. In the hyperbolic case, 
using the compactness argument gets complicated because of the lack of regularity over 
the whole time-space domain. 

Our approach is based on the fact that for a range of boundary operators, solutions 
improve smoothness dynamically, more precisely, they eventually become A:-times continu- 
ously differentiable for each particular k. We prove such kind of results in Section 2. Note 
that in some interesting cases the smoothing phenomenon was shown earlier in [12, 23]. 

This phenomenon allows us in Section 3 to work out a regularization procedure via 
construction of a parametrix. We here present a quite general approach to proving the 
Fredholmness for first-order dissipative hyperbolic PDEs and apply it to the periodic 
problems. Our Fredholm results cover non-strictly hyperbolic systems with discontinu- 
ous coefficients, but they are new even in the case of strict hyperbolicity and smooth 
coefficients. 

From a more general perspective, the smoothing effect and Fredholmness properties 
turn out to play an important role in the study of the Hopf bifurcation and periodic 
synchronizations in nonlinear hyperbolic PDEs [2, 21] via the Implicit Function Theorem 
and Lyapunov-Schmidt procedure [7, 14] and averaging procedure [6, 30]. 

From the practical point of view, our techniques cover the so-called traveling-wave 
models from laser dynamics [22, 28] (describing the appearance of self-pulsations of lasers 
and modulation of stationary laser states by time periodic electric pumping), population 
dynamics [9, 13, 34], and chemical kinetics [3, 4, 5, 36] (describing mass transition in 
terms of convective diffusion and chemical reaction and analysis of chemical processes in 
counterflow chemical reactors). 



2 Smoothing effect 

Here we describe classes of (initial-)boundary problems for first-order one-dimensional 
hyperbolic PDEs whose solutions improve their regularity in time. 
Set 

U T = {(x,t) : < x < 1, T < t < oo}. 

We address the problem 

(d t + a(x,t)d x + b(x,t))u = f(x,t), (2.1) 



u(x,0) = <p(x), (2.2) 
u J (0,t) = (Ru) J (t), l<j<m 

Uj(l,t) = (Ru)j(t), m < j <n 1 ' ' 

in the semi-strip n and the problem (2.1), (2.3) in the strip n_ 00 . Here u = (u±, . . . , u n ), 
f — • • • , fn), and ip = (ipi, . . . , ip n ) are vectors of real- valued functions, b = {bjk}^k=i 
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and a = diag(ai, . . . , a n ) are matrices of real- valued functions, and 1 < m < n are fixed 
integers. Furthermore, R is an operator mapping C (Ho) into C ([0, oo)) n , and similarly 
for R in Tl-oo- In Sections 2.1-2.3 we give examples of R as representatives of some classes 
of boundary operators ensuring smoothing solutions. 
In the domain under consideration we assume that 

a,j > for all j < m and aj < for all j > m, (2.4) 

inf | a,- 1 > for all j < n, (2.5) 

x,t 

and 

for all 1 < j ^ k < n, there exists pjk G C 1 such that bjk = Pjk{&j — cik)- (2.6) 

Note that all these conditions are not restrictive neither from the practical nor from the 
theoretical points of view. In particular, condition (2.4) is true in traveling-wave models 
of laser and population dynamics as well as chemical kinetics, where the functions Uj for 
j < m (respectively, m + 1 < j < n) describe "species" traveling to the right (respec- 
tively, to the left). Condition (2.5) means that all characteristics of the system (2.1) are 
bounded and the system (2.1) is, hence, non-degenerate. Finally, the condition (2.6) is 
a kind of Levy condition usually appearing to compensate non-strict hyperbolicity where 
the coefficients cij and for some j ^ k coincide at least at one point, say, (x , t ). In this 
case the lower-order terms with the coefficients bjk and bkj contribute to the system at 
(xo,to) longitudinally to characteristic directions (keeping responsibility for the propaga- 
tion of singularities), while in the strictly hyperbolic case we have a qualitatively different 
transverse contribution at that point. The purpose of (2.6) is to suppress propagation of 
singularities through the non-diagonal lower-order terms of (2.1). 

We will impose the following smoothness assumptions on the initial data: The entries 
of a, b, and / are C°°-smooth in all their arguments in the respective domains, while the 
entries of ip are assumed to be continuous functions only. 

Let us introduce the system resulting from (2.1)-(2.3) (resp., from (2.1), (2.3)) via 
integration along characteristic curves. For given j < n, x G [0, 1], and t G K, a charac- 
teristic of (2.1) is defined as the solution £ G [0, 1] >->■ u>j(^;x, t) G R of the initial value 
problem 

d^;x,t) = aj{ ^. Xit) y u Jj (x;x,t)=t. (2.7) 

Define 

Cj(f , x, t) = exp f f^A (rj, ujj(rj; x, t)) drj, dj(£, x, t) = ^^' ^'^ 

and 

Xj = if j < m and Xj — 1 if j > m. 
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Straightforward calculations show that a C 1 -map u : [0, 1] x [0, oo) — > R n is a solution to 
(2.1)-(2.3) if and only if it satisfies the following system of integral equations 



Uj(x,t) = (BSu)j(x,t) 




Jx i fc=i 



+ / dj{£,x,t)fj{£,Uj(£-,x,t))d£, j < n, 



(2.8) 



where 



(Bu)j(x, t) = Cj(Xj,X, t)Uj (Xj,UUj(Xj] x,t)) , 



(2.9) 



(Su)j(x,t) 



(Ru)j(t) if t > 0, 
(fj(x) if t — 0. 



(2.10) 



Here B is a shifting operator from <9n along characteristic curves of (2.1), while the 
operator S is used to denote the boundary operator on the whole dU . Similarly, a C 1 - 
map u : [0, 1] x R — > M n is a solution to (2.1), (2.3) if and only if it satisfies the system 
(2.8), where the definition of S is changed to S = R. 
This motivates the following definition: 

Definition 2.1 (1) A continuous function u is called a continuous solution to (2.1)-(2.3) 
in Ho if it satisfies (2.8) with S defined by (2.10). 

(2) A continuous function u is called a continuous solution to (2.1), (2.3) in n_oo if it 
satisfies (2.8) with S = R. 

Existence results for (continuous) solutions to the problems under consideration are 
obtained in [1, 15, 16, 18]. 

Definition 2.2 A solution u to the problem (2.1)-(2.3) or (2.1), (2.3) is called smoothing 
if, whatever k G N, there exists T > such that Uj G C k (ll T ) for all j < n. 

For the initial-boundary value problem (2.1)-(2.3) Definition 2.2 reflects a dynamic 
nature of the smoothing property stating that the regularity of solutions increases in time. 
The fact that the regularity cannot be uniform in the entire domain is a straightforward 
consequence of the propagation of singularities along characteristic curves. Moreover, 
switching from C k to C* +1 -regularity is jump-like; this phenomenon is usually observed 
in the situations when solutions of hyperbolic PDEs change their regularity (see e.g., [23, 



Note that, if the problem (2.1), (2.3) is subjected to periodic conditions in t, then 
Definition 2.2 implies that the smoothing solutions immediately meet the C°°-regularity 
in the entire domain. 



25, 27, 29]). 
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Definition 2.2 captures the general nature of the smoothing phenomenon for hyperbolic 
PDEs. A more precise information can be extracted from the proof of Theorems 2.3, 2.4, 
and 2.5 below: Reaching the C fc -regularity for solutions needs only a C fe+1 -regularity for 
a, b, and /. More exact regularity conditions for the boundary data, which also depend 
on k, can be derived from these proofs as well. These refinements are useful in some 
applications. 

Definition 2.2 can be strengthened by admitting worse regularities for the initial data. 
One extension of this kind, when the initial data are strongly singular distributions con- 
centrated at a finite number of points, can be found in [17]. In [17] we used a delta-wave 
solution concept. Another result in this direction [19, 20] concerns periodic problems and 
uses a variational setting of the problem (see also Theorem 3.2 (ii)). In [19, 20] we get an 
improvement of the solution regularity from being functionals to being functions. 

In what follows we demonstrate the smoothing effect on generic examples of large 
classes of boundary operators and show which kinds of problems can be covered by our 
techniques. Our approach to establishing smoothing results is based on the consideration 
of the integral representation of the problems and observation that the boundary and the 
integral parts of this representation have different influence on the regularity of solutions. 
Our main idea is to show that the integral part has a self-improvement property, while 
the boundary part can in many interesting cases be not responsible for propagation of 
singularities. The latter contrasts to the case of the Cauchy problem where the solutions 
cannot be smoothing as the boundary term all the time "remembers" the regularity of 
the initial data. It is worthy to note that in the case of the problem (2.1)-(2.3) in n the 
domain of influence of the initial conditions is determined by both parts of the integral 
system and is in general infinite. This makes the smoothing effect non-obvious. 

2.1 Classical boundary conditions 

Here we specify conditions (2.3) to 

Uj(0,t) = hj(t), l<j<m, , 2n s 
Uj(l,t) = hj(t), m < j < n. 

and consider the problem (2.1), (2.2), (2.11). 

Theorem 2.3 Assume that the data aj, bjk, fj, and hj are smooth in all their argu- 
ments and Lpj are continuous functions. Assume also (2.4), (2.5), and (2.6). Then any 
continuous solution to the problem (2.1), (2.2), (2.11) is smoothing. 

Note that in the case of smooth classical boundary conditions (2.11), the domain of 
influence of the initial data (p(x) on U{ for every i < n in general is unbounded (due to 
the lower-order terms in (2.1)). In spite of this, the influence of the initial data on the 
regularity of u becomes weaker and weaker in time causing the smoothing effect. 

Proof. Suppose that u is a continuous solution to the problem (2.1)-(2.3) and show 
that the operator of the problem improves the regularity of u in time. The idea of the 
proof is similar to [17]. 
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We start with an operator representation of u. To this end, introduce linear bounded 
operators D, F : C (TT )" ->■ C (n )" by 

(Du)^ (x, t) = - f d j (£,x,t)^b jk (£,uj j (£;x,t))u k (t,uj j (£;x,t))dt, 

Jxj k = l 

Note that Ff is a smooth function in x,t. In this notation the integral system (2.8) can 
be written as 

u = BSu + Du + Ff. (2.12) 

It follows that 

u = BSu + (DBS + D 2 )u + (/ + D)Ff. (2.13) 

In the first step we prove that the right hand side of (2.13) restricted to IT^ for some 
Ti > is continuously differentiable in t. The C 1 (iTrJ -regularity of u will then follow 
from the fact that u satisfies (2.1) in the distributional sense. By the assumption (2.5), 
we can fix a large enough 7\ > such that the operator S in the right-hand side of (2.13) 
restricted to IT^ does not depend on ip and, hence, Su = Ru = h, where h = (hi, . . . , h n ). 
We therefore arrive at the equality 

u | n = Bh + DBh + D 2 u + (I + D)Ff, (2.14) 

where u\ji t denotes the restriction of u to IT^. By the regularity assumption on a, b, f, 
and h, the function Bh + DBh + D 2 u + (I + D)F f is smooth. We are reduced to show 
that the operator D 2 is smoothing, more specifically, that D 2 u is C 1 -smooth in t on IT^. 
Fix a sequence u l E C 1 (n ) such that 

u l ->■ u in C (n ) n as I ->■ oo. (2.15) 

Under convergence in C (fi) n here and below we mean the uniform convergence on any 
compact subset of f2. Then D 2 u l — > D 2 u in C (n ) as well. It suffices to proof that 
d t [D 2 u l ] converges in C(n Tl ) n as/ — > oo. Given j < n, consider the following expression 
for (D 2 u l ) . (x, t), obtained after the change of the order of integration: 



(D 2 u l ).(x,t) (2.16) 

n n „ x „ x 

= X X / / d 3ki(€, V, x, t)b jk (£, Uj(£; x, t))u\(r], u k (r]; £, Ujfc x, t)))d£dr} 

t-i ,-i JO J n 



k^j i^k 



with 

djkM, V, x, t) = dj(£, x, t)d k (r], f , x, t))b ki (r], u k (r]\ f , x, t))). 
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It follows that 

d t [(D 2 u l ).(x,t) 

n n 

= EE 



fc=l i=l 



+ EX 

fc=l i=l 



Jr) 



dt [djkiit, V, x, t)b jk (£, x, t))] ufa, u k (r]; f , Uj(£; x, t)))d£dr) 



Jr) 



djkiH, V, x, t)b jk (€, x, t)) 



x d 3 u k (r); f , ujj(£; x, t))dtWj(& x, t)d 2 u\{rj, uj k (r); f , u^(£; x, t)))d£dr), 



(2.17) 



where d r g here and below denotes the derivative of g with respect to the r-th argument. 
Our task is therefore reduced to show the uniform convergence of all integrals in the 
second summand, whenever (x,t) varies on a compact subset of IL^. For this purpose we 
will transform the integrals as follows. Using the assumption (2.6) and the formulas 



d x ujj(£;x,t) 



1 



exp 



• r 'ha. 



a,j(x, t) 

dtWj(£;x,t) = exp J (^^j x,t))dr), 



(v,Uj(v,x,t))dr), 



(2.18) 
(2.19) 



we get 



Jr) 



djki(€, V, x, t)b jk (£, ujj(£; x, t)) 
x d 3 u k (ri; f , Wj(Z; x, t))d t Uj(£; x, t)d 2 u\(r], u k (r); £, Wj(£; x, t)))d£dr] 
djki(t V, x, t)d 3 u k (ir, f , x, t))dtUj(Z; x, t) 



Jr) 



PX PX 

= d jki (^,r],x,t)dtuj j (^x,t)(a k a j )(^,uj j (^x,t)) 

JO Jr) 

bjk ) (£, ; x, t)) (d ( uty (77, uj k (r}; f , Ujfc x, t)))d£dr) 

a k — a j J 

X fX 

djki(€, V, x, t) (a c u{) (?/, w fc (7/; f , ; x, t)))d£dr] 
d^djki(C, V, x, t)u\ (77, u k (r); f , Uj(£; x, t))) d^drj 



Jr) 



+ 



f 

JO 



Jr) 



djki{£, V, x, t)u\ (77, a; fc (r/; f , c^-(£; x, *))) dr/. 

J £=?} 



Here 



»7, x, t) = d jki (£, rj, x, t)dtWj(Z; x, t) 



Now, the desired convergence follows from (2.15). 

In the second step we prove that there exists T 2 > T\ such that d t u restricted to 
Ut 2 is C 1 -smooth in t on Ht 2 - Once this is done, we differentiate (2.1) in t and get 
d 2 t u G C (n T2 ) ; differentiating (2.1) in x, we get d^u G C (n T2 ) . We will be able 
to conclude that u G C 2 (Rt 2 ) , as desired. To prove the existence of T 2 , let v = d t u. 
Differentiation of (2.1) formally in t leads to 

n n 

(d t + a j d x )v j + ^ b JkVk + ^2 d t b JkUk + d t ajd x Uj = d t fj. 

k=l k=l 

Combining this with (2.1), we obtain 

(d t + a.jd x )vj + ^2 h ikVk —Vj 

k=i a i 

n q / n \ 

= d tfj ~ Yl d ^ u k + — Yl h i kU i ~~ fi = G i(-f'v d tfj,u). (2.20) 
k=i a i \k=i J 

Here, for each j < n, Gj is a linear function with smooth coefficients. Set 

. H fbjj d t cij\ , , , ~ . . . Cj(£,x,t) 

Cj(Z, x,t)=exp - — / (77, Uj(ri; x, t)) drj, dj{£, x, t) = J 

Jx \ a j a j J fl ils) w jlsi x i l )) 

and introduce three linear operators B,D,F : C (Ho) — > C (Ho) by 

{Bu^J ( x : t) = Cj(Xj,X, t)Uj (Xj,Uj(Xj] X, t)) , 

[Du] (x,t) = - djfox^^bjkfaujfcx^Ukfaujfcx^dt, 
(Ff),(x,t) = f\l^.xJ)J)^.^:xJ)Uli. 

J Xj 

Similarly to the above, our starting point is that there is T 2 > T\ such that v satisfies the 
following operator equation resulting from (2.20): 

V |n T2 = Bti + Dv + FG(f,d t f,u), 

and, hence, the equation 

v|n T2 = Bti + DBti + D 2 v + (/ + D)FG(f, d t f, u), (2.21) 

where G = (G 1: . . . , G n ) and h! = (h[, . . . , h' n ). Again, due to the assumption (2.5), we 
can fix T 2 > Ti such that the right-hand side of (2.21) does not depend on u and v in 
IT \ Htx- Due to Step 1, the function (I + D)FG(f, dtf, u) then meets the -regularity. 
We are thus left to show that the operator D 2 is smoothing in the above sense. As D is 
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exactly the operator D with Cj and dj replaced by the smooth functions 5j and dj, the 
desired smoothing property of D 2 follows from the proof of the smoothness of D 2 and the 
fact that D 2 v in (2.21) does not depend on v in II \ IT^. 

Proceeding further by induction, assume that, given r > 2, there is T r > such 
that u G C r (HT r ) and prove that u meets the C r+ ^regularity in t on Hr r+1 for some 
T r+ i > T r . Set w = d r t u and write our starting operator equation for w resulting from 
(2.1), (2.3) after r-times differentiation in t: 

iy| n = Bh^ + DBh^ + D 2 w + (I + D)FG(f, d t f, d r t f, u, D 1 u, D^u) 

(2.22) 

with G being a vector of linear functions with smooth coefficients. Similarly to the above, 
fix T r+ \ > T r such that the right-hand side of (2.22) does not depend on u, D 1 u, . . . , D r ~ l u, 
and w in IT\ITt,,. This ensures that the last two summands in (2.22) are C^-functions. The 
first two summands are C^-smooth by the regularity assumptions on the data. Finally, 
the C r+1 (ll Tr+1 ) -regularity of u follows from the previous steps of the proof and suitable 
differentiations of the system (2.1). □ 

Theorem 2.3 can be extended over the boundary operators of the following kind (both 
linear and nonlinear). Given T > 0, in the domain T1 T let us consider the problem (2.1)- 
(2.3) with bij = for all i ^ j (i.e., the system (2.1) is decoupled) and with (2.2) replaced 
by u(x,T) = (p(x) (the initial values are given at t = T). This entails that the domain 
of influence of <p depends only on the boundary conditions. For the latter it is supposed 
that, whatsoever T > and <p(x), the function tp(x) has a bounded domain of influence on 
u. In other words, for any decoupled system (2.1), if ip(x) has a singularity at some point 
x G [0, 1], then this singularity expands along a finite number of characteristic curves, 
and this number is bounded from above uniformly in i £ [0,1]. This class of boundary 
operators is in detail described in [17], where the necessary and sufficient conditions 
for smoothing solutions are given. The results of [17] generalize the smoothing results 
obtained in [12, 23] for the system (2.1) with time-independent coefficients and (a kind 
of) Dirichlet boundary conditions. 

2.2 Integral boundary conditions in age structured population 
models 

Here we address another class of boundary operators admitting smoothing solutions. 
Though it covers a range of (partial) integral operators, we illustrate our smoothing 
result with an example from population dynamics. 

Integral boundary conditions are usually used in continuous age structured population 
models to describe a fertility of populations. Let u(x, t) denotes the density of a population 
of age x at time t. Then the dynamics of u can be described by the following model (see, 
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e.g. [10, 24, 34] and references therein): 



(d t + d x + n) u = o, (x, t) e n , (2.23) 

u(x,0) = (p(x), x e [0,1], (2.24) 
u(0,t) = /iQf j(x)u(x,t) dxj , teR, (2.25) 

where > is the mortality rate of the population and the functions h and 7 describe the 
fertility of the population. Not losing potential applicability to the topic of population 
dynamics, h and 7 are supposed to be C°°-smooth functions. The integral in (2.25) is 
a kind of the so-called "partial" integral, since u depends not only on the variable of 
integration x, but also on the free variable t. By this reason the right-hand side of (2.25) 
is not smoothing. Nevertheless, it turns out that it is regular enough to contribute into 
smoothing solutions. 

Theorem 2.4 Assume that h and 7 are C 00 -smooth functions and ip is a continuous 
function. Then any continuous solution to the problem (2.23)-(2.25) is smoothing. 

Proof. It suffices to show the smoothing property starting from large enough t. 
Therefore, we can use the notation: 



(Ru)(t) = h (^j r ){x)u{x,t) dx 



u(£;x,t) = t + ^-x 
c(£,x,t) = c(Z,x,t) = 
(Bu)(x,t) = (Bu)(x,t) = e-^«(0,t-x), 

the latter two being introduced for all large enough t. Integration along the character- 
istic curves implies that any continuous solution to (2.23)-(2.25) satisfies the operator 
equations u = BRu and u = Bu and, hence, 

u = BRBu (2.26) 

whenever t > 7\, where Ti is chosen to be so large that the operator BRB moves away 
from the initial boundary (the right-hand side of (2.26) does not depend on ip). Since 

(BRBu)(t) = e-^h(^j\(Qe-*u(0,t-x-OdS) 

c h QT -y(t-x- T)e^ x - t+T) u{<d, r) dr^j , 



we obtain the C^-smoothness of BRBu and, hence, of u on n Tl . The C 1 -smoothness of 
u on n Tl now follows from (2.23). 
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Proceeding similarly to the proof of Theorem 2.3, in the second step we consider the 
following operator equation with respect to v — d t u, obtained after differentiation of (2.23) 
and (2.25) in t and integration along characteristic curves: 



ln T , 



= Bd+RBv, 



(2.27) 



where 



(d t Rv)(t) — h! ( / r y(x)u(x : t) dx ) / r y(x)v(x, t) dx 



and T 2 > Ti is fixed to satisfy the property that the right-hand side of (2.27) does not 
depend on u and v in n \ IT Tl . It follows that 



ln T , 



= e-^ti 



f 7(fl«(& t - x) d£\ £ -y(Z)e->*v(0, t - x - d£ 

l(Ou(£,t-x)dA [ 1 (t-x-T)e^ x - t+r ^v(0,r)dT. 



t-x-l 



To conclude that v G C\ (n^)™, it remains to note that u under the first integral in the 
right-hand side meets the C^-regularity, while the second integral gives us the desired 
smoothing property. 

In general, given T r for r > 2, we choose T r+ i > T r by the argument as above and for 
w = &[u have the equation 



H n „ =Bd r t RBw, 

I T' + 1 



(2.28) 



where 



(%Rw){t) = h'^ 



r y(x)u(x : t) dx J / r y(x)w(x, t) dx 



+ 
+ 



d 



r-l 



dr- 1 
d r ~ 2 



dt 



r-2 



ti 
ti 



j(x)u(x, t) dx 



j(x)d t u(x, t) dx 



r y(x)u(x, t) dx J / r y(x)d t u(x : t) dx 



Substituting the latter into (2.28) and changing variables under the integral of w similarly 
to the first two steps, we get the desired smoothing property for w. This completes the 
proof. □ 



2.3 Dissipative boundary conditions 

Now we switch to boundary conditions having dissipative nature and fitting the smoothing 
property. A large class of dissipative boundary conditions for hyperbolic PDEs is described 
in [8]. 
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To give an idea of the smoothing effect in this case, consider the following specification 
of (2.1): 

Uj{0,t) = hj(z(t)), l<j<m, 
Uj(l,t) = hj(z(t)), m < j < n, 

with 

z(t) = («i(M), . . . ,u m (l,t),u m+1 (0,t), . . .,u n (0,t)) . 

In the domain n_oo we address the problem (2.1), (2.29) subjected to periodic boundary 
conditions 

u(x,t + 2n) = u(x,t). (2.30) 

The problems of this kind appear in laser dynamics and chemical kinetics (in Section 3 
we investigate a traveling-wave model of kind (2.1), (2.29), (2.30) from laser dynamics). 
Within this section, using the standard notation for the (sub) spaces of continuous func- 
tions, we mean that the functions have additional property of 27r-periodicity in t. Write 

n 

h'(z)=V z h J (z)> IWI=EM(*)|, = {W*) }l k =v 

k=l 

Theorem 2.5 Assume that cii, bij, fa, and hi are smooth functions in all their arguments 
and the conditions (2.4)-(2.6) are fulfilled. Moreover, the functions ai, b^, fa are supposed 
to be 2n -periodic in t. If 

|^(z)| expj^ 3 (^j foc^faM))^} < 1 (2.31) 

and 

|/»;(z)|oq> jjT (^-^) (V^j(v,x,t))dr^<l (2.32) 

for all j < n, x e [0, 1], t e K, and z e ffi™, then any continuous solution to the problem 
(2.1), (2.29), (2.30) is smoothing. 

Proof. Any continuous solution to the problem (2.1), (2.29), (2.30) in n_oo fulfills 
(2.12) with S = R and also satisfies the equation 

u = Bu + Du + Ff (2.33) 

where the boundary conditions are not specified. Substituting (2.33) into (2.12), we obtain 

u = BRu+(DB + D 2 )u+(I + D)Ff. (2.34) 

We first show the bijectivity of I — BR e C [C] (n^)"). On the account of (2.19) and 
the definition of B given by (2.9), we have 

(BRu)j(x,t) = Cj(xj,x,t)hj (z(u>j(xj;x,t))) = Cj(xj,x,t)hj(0) 

+ exp|y ^— ^ (rj, ujj(r); x, t)) drj^ J h'j (az(cOj(xj;x,t)) da ■ z(ujj(xj;x,t)) 
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and 

d t [(BRu)j(x, t)} = (d t Cj) (xj,x, t)hj (z(u j (x j ; x, t))) 

+ h 'j ( z (uj(xj] x, t)) ■ z'iujjixj] x, t)) exp | J - ^J-^ (rj, Uj(r)] x, t)) drj 

where • denotes the scalar product in IR n . Taking into account (2.30), the bijectivity of 
I— BR E C (p] (n_oo) ) now follows from the contractibility conditions (2.31) and (2.32). 

Now we claim that the operators DB and D 2 in (2.34) are smoothing. The latter is 
smoothing by the proof in Theorem 2.3. Similar argument works also for DB. Indeed, 
by the definition of the operators D and B we have 

(DBu l ).(x,t) (2.35) 



= z2 d i^> x ' > x > l )) c k{xk, £, x, t))u[(x k ,uj k (x k ; f , x, t)))d£, 

k=l 

where the sequence u l is fixed to satisfy (2.15) with n replaced by n_oo. To show that 
dt [DB-u'j converges uniformly on n_oo, we transform the integrals in (2.35) like to the 
case of D 2 , that is, we differentiate (2.35) in t, use (2.6), and integrate by parts. In this 
way we get the smoothing property for DB. Turning back to the formula (2.34) and using 
in addition the fact that (I + D)Ff is C°°-smooth, we can rewrite (2.34) in the equivalent 
form 

u = (I - BR)- 1 [(DB + D 2 )u + (I + D)Ff] , 

thereby reaching the C\ -regularity for u. Afterwards, the C 1 -regularity of u is a straight- 
forward consequence of the system (2.1). 

Proceeding similarly to the proof of Theorem 2.3, we come to the formula for v = d t u: 

v = (I- BR',)- 1 [(DB + D 2 )v + (/ + D)FG(f, dj, uj\ , 

where R' z y = h'(z)y. It follows from the bijectivity of / — BR' Z e C (C] (n_oo) n ), which 
we have by the same conditions (2.31) and (2.32), and the C^-regularity of DB + D 2 and 
(I + D)FG(f,d t f,u). Therefore, v E C\ (IL^)", what entails u E C 2 (n^)". It follows 
by (2.1) that uEC 2 (TT^)" '. 

To complete the proof, we proceed by induction on the order of regularity of u. Assume 
that u E C r (iLoo) for some r > 2 and prove that u E C r+1 (n_ 0O ) . Our starting 
formula for w = &[u is as follows: 

w = (I — BR'^ 1 \(DB + D 2 )w 

+ (/ + D)FG(f, d t f, d r J, u, d t u, dr\) 

+ m- x R' z d t z + m t - 2 (R' z d t z) 

where d r t ~ x R! z = {dl~ 1 (dkhj(z))} n . By the regularity assumptions on the data and the 
induction assumption, the last three summands in the square brackets are C^-functions. 
Using in addition our smoothing argument for DB + D 2 and the regularity properties of 
(I — BR^)- 1 , we arrive at the desired conclusion. □ 
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3 Fredholm solvability of periodic problems 

In [19, 20] we suggested an approach to establishing the Fredholm property for first- 
order hyperbolic operators. This is done by construction an equivalent regularization 
in the form of a parametrix. The construction is, implicitly but essentially, based on 
the smoothing effect investigated in Section 2. Consider the first-order one-dimensional 
hyperbolic system 

(d t + a{x)d x + b{x))u = f(x,t), x e (0,1), (3.36) 
subjected to periodic conditions (2.30) and reflection boundary conditions 

n 

Uj (p,t)= l<j<m, 

k=m+l 

m 

U A 1 ^) = Y r jk u k(l,t), m<j<n. 
k=i 

Here r® k and r l - k are real numbers and the right-hand sides fj : [0, 1] x K. — > R are supposed 
to be 27r-periodic with respect to t. 

The main result of this section states that the system (3.36), (2.30), (3.37) is solvable 
if and only if the right hand side is orthogonal to all solutions to the corresponding 
homogeneous adjoint system 

— d t u — d x (a(x)u) + b T (x)u = 0, x G (0, 1), 

subjected to periodic conditions (2.30) and adjoint boundary conditions 

m 

a,-(0)?ij(0, t) = -^r° kj a k (0)u k (0,t), m < j < n, 

k=l 

n 

aj(l)uj(l,£) = - Y r lj a k(l)uk(l,t), 1 < J < m 

k=m+l 

We will present our result in three steps. First we introduce appropriate function 
spaces for solutions. Then we decompose the operator of the problem into two parts, only 
one being responsible for propagation of singularities. Finally, based on this decomposition 
and the smoothing property, we construct a parametrix thereby establishing the Fredholm 
solvability. 

Choosing the function spaces, note that the problem (3.36), (2.30), (3.37) describes 
the so-called traveling- wave models from laser dynamics [22, 28]. From the physical point 
of view, it is desirable to allow discontinuities in the coefficients and the right hand side 
of (3.36). This entails that the spaces of solutions should not be too small. On the 
other hand, they should not be too large, in order to admit embeddings into an algebra of 
functions with pointwise multiplication of its elements. The last property is important for 
potential applicability of our results to nonlinear problems, like describing such dynamic 
phenomena as Hopf bifurcation and periodic synchronizations. Finally, the solution spaces 
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(3.37) 



(3.38) 



capable to capture the Fredholm solvability need to have optimal regularity with respect 
to the function spaces of the right-hand side. 

We now describe the scale of spaces V 1 (for the solutions) and W 1 (for the right-hand 
side) meeting all these properties. For 7 > 0, let W 1 denote the vector space of all locally 
integrable functions / : [0, 1] x R — > W 1 such that f(x,t) = f (x, t + 2n) for almost all 
x G (0, 1) and t G R and that 



2 

W 1 



Ed 



2?r 



f(x,t)e- tst dt 



dx < 00. 



(3.39) 



Here and in what follows || ■ || is the Hermitian norm in C n . It is well known that W 1 is 
a Banach space with the norm (3.39); see, e.g. [11], [31, Chapter 5.10], and [33, Chapter 
2.4]. 

Furthermore, for 7 > 1 and a G L°° ((0, l);M n ), where M„ denotes the space of real 
n x n matrices, with ess inf \a,j\ > for all j < n we will work with the function spaces 

W = | M g W 1 : d x u G W°, d t u + ad x u G W 7 } 

endowed with the norms 



u 



u 



|2 



+ \\d t u + ad x u\ 



Remark that the space U 1 depends on a and is larger than the space of all u G W 1 such 
that dtu G W 1 and d x u G W 1 (which does not depend on a). For u G U 1 there exist 
traces u(0, •), u(l, •) G L 2 oc (R;R n ) (see [20]), and, hence, it makes sense to consider the 
closed subspaces in U y 



V = {ueW : (3.37) is fulfilled}, 
V~i = { u G W : (3.38) is fulfilled}. 

Our next task is to decompose the operator of our problem into two parts in order 
to single out the part, denoted below by A, which is bijective and at the same time is 
responsible for the propagation of singularities. If this decomposition is optimal, then 
after a regularization procedure the other part becomes smoothing and therefore meets 
the compactness property. Let 

b° = diag(6n, b 22 , ■ ■ ■ , b nn ) and b L = b - b° 

denote the diagonal and the off-diagonal parts of the coefficient matrix b, respectively. 
Let us introduce operators A G C(V^; W"'), A G £(F 7 ; W~<), and B, B G C{W~*) by 

Au = d t u + ad x u + b°u, 

Au = —d t u — d x (au) + b°u, 

Bu = 

Bu = {bYu. 
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Remark that the operators A, B, and B are well-defined for aj, bj k G L°°(0, 1), while A is 
well-defined under additional regularity assumptions with respect to the coefficients aj, 
for example, for a,- G C 0,1 ([0, 1]). Note that the operator equation 

Au + Bu = f 

is an abstract representation of the periodic-Dirichlet problem (3.36), (2.30), (3.37). 
Finally, for s G Z we introduce the complex (n — m) x {n — m) matrices 



R s 

where 



'jl'lk 



j,k=m+l 



aAx)= L W) dv ' m= L b ^M dy - 

The following theorem states, first, that the pair of spaces (V 1 , W 1 ) gives an optimal 
regularity trade-off between the spaces of solutions and right-hand sides and, second, 
that A meets the bijectivity property. The second desirable property for A of being an 
optimal operator responsible for propagation of singularities will be a consequence of our 
Fredholmness result. 

Theorem 3.1 [20] For every c > there exists C > such that the following is true: If 
aj, bjj G L°°(0, 1) and essinf \aj\ > c for all j — 1, ... , n, (3.40) 

n m n n m 

jr'=l j=l k=m+l j=m+l k=l 

and 

| det(7 - R s (a, b°, r))\ > c for all s G Z, (3.41) 
then for all 7 > 1 the operator A(a, b°) is an isomorphism from V 1 {a, r) onto W 1 and 

\\^{.tt,^y 1 \\ C {W-i-yi(a > r)) < C. 

Let 

(J,«) L 2 = — y J (f(x,t),u(x,t)) dxdt 

denote the scalar product in the Hilbert space L? ((0, 1) x (0,27r);M n ) and (•, •) denote 
the Euclidean scalar product in M n . We are prepared to formulate the main result of this 
section. 

Theorem 3.2 [20] Suppose that conditions (3.40) and (3.41) are fulfilled for some c > 0. 
Suppose also that 

for all j 7^ k there is Cjk G -6^(0, 1) such that , , 

a k (x)b jk (x)a = c jk (x)(aj(x) - a k (x)) for a. a. x G [0, 1]. 
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Then the following is true: 

(i) The operator A + B is a Fredholm operator with index zero from V 1 into W 1 for 
all 7 > 1, and kei(A + B) = {u e V 1 : (A + B) u = 0} does not depend on 7. 

(ii) (smoothing effect) If a e C ' 1 ([0, 1]; M n ) ; then ker(A + B)* = ker(A + B) and 

{{A+B)u : u e V 1 } = [/ e W~< : (/,m) L 2 = /or a// w G ker(„4 + £)} , 

where ker(^l + = V 7 : (^4. + = 0} does not depend on 7. 

Theorem 3.2 (ii) states that the kernel of the adjoint operator is actually defined 
on the classical function spaces. In other words, the kernel has much better regularity 
than ensured just by the formal definition of the adjoint operator. Here we encounter a 
smoothing effect for the solutions (of the adjoint hyperbolic problem), that are originally 
functionals. The proof of this effect in [19, 20] uses completely different techniques, based 
on a functional-analytic approach. 

Finally, we outline the proof of Theorem 3.2 (i). As mentioned above, we construct a 
parametrix to the operator of the problem. By Theorem 3.1, the zero-order Fredholmness 
of the operator A + B £ £(y 7 ; W 1 ) is equivalent to the zero-order Fredholmness of the 
operator I + BA^ 1 G CiW 1 ). Furthermore, we use the following Fredholmness criterion 
(see also [32, Theorem 5.5] or [35, Proposition 5.7.1]). 

Lemma 3.3 [19] Let I denote the identity in a Banach space W. Suppose that V e C{W) 
and V 2 is compact. Then I + V is Fredholm. 

Setting V = BA' 1 E £(H/ 7 ), we prove that V 2 e C{W~<) is compact (while V alone 
can hardly be compact, being a type of a partial integral operator). This actually means 
that T> 2 has smoothing property. In fact, T> 2 is basically the same as the operator D 2 , 
that we used in the proof of Theorem 2.3. 

Since / - V 2 = (I - D)(I + D) = (I + D)(I - D), the operator / - V is a parametrix 
of / + V. It follows that the operator A + B admits an equivalent regularization in the 
form of the right parametrix A~ l {I — BA~ V ). 
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